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Abstract 

In this review we present recent theoretical results concerning investigations of single subsurface 
defects by means of a scanning tunneling microscope (STM). These investigations are based on 
the effect of quantum interference between the electron partial waves that are directly transmitted 
through the contact and the partial waves scattered by the defect. In particular, we have shown 
the possibility imaging the defect position below a metal surface by means of STM. Different types 
of subsurface defects have been discussed: point-like magnetic and non-magnetic defects, magnetic 
clusters in a nonmagnetic host metal, and non-magnetic defects in a s-wave superconductor. The 
effect of Fermi surface anisotropy has been analyzed. Also, results of investigations of the effect 
of a strong magnetic field to the STM conductance of a tunnel point contact in the presence of a 
single defect has been presented. 

PACS numbers: 61. 72. J Point defects and the defect clusters, 73.40.Cg Contact resistance, contact potential, 
73.63.Rt Nanoscale contact, 74.50.+V Tunneling phenomena: single particle tunneling and STM, 73.23.-b 
Electronic transport in mesoscopic systems, 72. 10. Fk Scattering by point defects, dislocations, surfaces, and 
other imperfections (including Kondo effect) 
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I. INTRODUCTION 

About three decades following its invention [l[, scanning tunnelling microscopy (STM) 
has proved to be a superbly valuable tool for investigating surfaces on the atomic scale. Along 
with a mapping of the conductor's surface, the STM enables observing many phenomena, 
among which electron scattering by single surface defects (impurity atoms, adatoms, or 
step edges). There are hundreds of papers that are devoted to investigations of surface 
defects by STM. In this paper we do not aim to review all of them and confine ourselves 
to briefly mentioning the main directions of researches in this field. Our attention will be 
mainly focused on interference effects in STM conductance caused by defects sitting below 
the surface. 

Electron scattering by defects leads to quantum-interference patterns in the local electron 
density of states around the defects (Friedel oscillations (2}). For more than thirty years 
Friedel oscillations have remained a theoretical prediction that could be seen only in theory 
textbooks j3[ . The appearance of the STM has enabled the visualization of these oscillations, 
which manifest themselves as oscillations of the differential tunneling conductance, G = 
dl/dV, around defects on the surface. 



2 



First standing wavelike patterns in the STM conductance in the vicinity of defects were 
observed by Crommie et al. [4] on a Cu(lll) surface and by Hasegava et al. [5| on a 
Au(lll) surface. At the (111) surface of the noble metals Cu, Ag, and Au the electrons of 
the surface states form a quasi-two-dimensional nearly-free electron gas having an isotropic 
dispersion law 0]. When scattered from step edges or adatoms the surface states form stand- 
ing waves which result in an oscillatory dependence of the tunneling conductance measured 
as a function of the distance between the STM tip and the defect, r . The period of the 
conductance oscillations Atq = 27r/2/c| D is set by twice the Fermi wave vector, 2kp D (kfr is 
a two-dimensional vector in the plane of the surface). 

The circular 2D Fermi contour of the electrons at the (111) surface of noble metals results 
from the fact that the layer of surface atoms actually corresponds to one of the close-packed 
stackings on which the face-centered cubic structure is based. Generally, for less closely 
packed surfaces and conductors having a complicated crystallographic structure a 2D Fermi 
contour is anisotropic, i.e. the absolute value of the vector kjP depends on its direction. 
The Fourier transform (FT) of the standing wave pattern provides an image of the Fermi 
contour. Anisotropic Friedel-like oscillations have been observed by FT-STM on Cu(llO) 
surfaces [3], Be Js}, and ErSi2 Particularly, in Ref. j3] the contour related to the 'neck' 
of the bulk Fermi surface for Cu (110) surface has been imaged. 

Magnetic adatoms on non-magnetic host metal surfaces are of special interest as they 
produce a characteristic many-body resonance structure in the differential conductance near 
zero voltage bias attributed to the Kondo effect [IMi. The shape of the resonance in the 



differential conductance is usually asymmetric and is described by a Fano line shape [14M16 
The surface electron waves carry information on the magnetic impurity and by focussing 
the waves it has been possible to create a mirage image of the impurity [17J (for review, see 
[l8|). The interesting phenomenon of an orbital Kondo resonance was observed by STM in 
Ref. [19]. It was found that STM images of the Cr(001) surface show cross-like depressions 
centered around the impurities corresponding to the orbital symmetry of two degenerate 
d xz , dyz surface states 19l . 



The investigation of defects near the surface of unconventional superconductors by STM 
is a way to determine the symmetry of the order parameter. The effect of single Zn defects 
on the superconductivity in high-T c superconductors was investigated in Ref. [2(j, and the 
manifestation of d-wave pairing symmetry was observed in the quasibound state near the 
defect. In Ref. [2lJ a bound state near a magnetic Mn adatom on the surface of supercon- 
ducting Nb was observed by STM. 

An effective way to enhance the STM sensitivity to such oscillation effects is to use a 
superconducting tip (22|. In Ref. 23] it was demonstrated that the amplitude of conductance 
oscillations is significantly enhanced when a superconducting tip is used, and when the 
applied bias is close to the gap energy of the superconductor. 

The applicability of STM can be extended to the study of magnetic objects on the surface 
of a conductor when a magnetic material is used for the STM tip such that the electric 
current is spin polarized (SP) (for review of SP-STM see j24j). For example, the precession 
of a magnetic moment of clusters of organic molecules on a surface gives rise to a time 



modulation of the SP-STM current, from which the g - factor can be found |25l . [26|. The 
possibility to probe magnetic properties of nanostructures buried beneath a metallic surface 
by means of local probe techniques is discussed in Ref. [27| . It has been shown that those 
properties can be deduced from the spin-resolved local density of states above the surface 
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STM spectroscopy also provides access to information on the structure of the metal 
below the surface in both semiconductors and metals. Crampin [28] proposed to utilize the 
surface states for imaging subsurface impurities. However, the exponential decay of the 
wave function amplitude into the bulk limits the effective range to the topmost layers only 
and bulk states form a good alternative for detecting defect positions. The principle of 
imaging subsurface defects is based on the influence on the conductance caused by quantum 
interference of electron waves that are scattered by defects and reflected back by the contact. 
This effect was explored for investigating subsurface Ar bubbles submerged in Al [29] and Cu 
Hj, and Si(lll) step edges buried under a thin film of Pb 31|. In these experiments, bulk 



electrons are found to be confined in a vertical quantum well between the surface and the 
top plane of the object of interest. The observation of interference patterns due to electron 
scattering by Co impurities in the interior of a Cu sample was reported Refs. [3J 




Reviews of STM theory can be found in Refs. [34], |35|. The papers listed in [34J, [35|, in 
which the conductance of a tunnel contact of small size has been analyzed theoretically, must 
be complemented by reference to the fundamental paper of Kulik, Mitsai and Omelyanchouk 
[36| published in 1974. In this paper the authors obtained, on the basis of rigorous quantum- 
mechanical considerations, an analytical formula for the conductance of a junction between 
two metal half-spaces separated by an inhomogeneous tunnel barrier of low transparency. 
Their result is valid for arbitrary values of the applied bias and for arbitrary dependence 
of the tunnelling probability on the coordinates in the plane of the interface between the 
metals. As a special case, the general formula for the contact resistance can be applied to an 
inhomogeneous tunnel contacts having a characteristic diameter smaller than electron wave 
length, which is suitable to describe STM conductance. Recently, electron tunnelling through 
a randomly inhomogeneous barrier of arbitrary amplitude has been analyzed theoretically 



in Refs. 37, 38 



The theoretical descriptions of STM conductance oscillations due to electron scattering 
by single defects in the majority of papers is based on the assumption that the tunnelling 
conductance measured by the STM tip is proportional to the local density of states (LDOS) 
u(r) of the sample (see, for example, [l8, 28, 39, 4o||) as for a planar tunnel junction j4l[. For 
the scattering of electron surface states this assumption is quite reasonable, but for electron 
scattering in the bulk of the sample it can not be used. The LDOS in the vicinity of defects 
in the bulk is critically modified by electron reflections off the surface of the conductor, at 
r G S, and differs from Friedel oscillations of the LDOS in an infinite conductor with a single 
scatterer [3|. In the limit of zero tunnelling probability we have u(r 6 S) = 0. Further, the 
conductance oscillations are formed only by "tagged" electrons, which tunnel through the 
contact and are scattered back by the defect, while a "halo" of Friedel oscillations around 
the defect is due to all scattered electrons. In general, there are no other periods in the 
interference effects but the period of Friedel oscillations Ar = 27t/2/cf (kp is a Fermi wave 
vector) and the analysis in Ref. [33j of the experimental data in terms of a bulk LDOS 
seems to be qualitatively correct (42|. However, the calculation of amplitudes and phases 
of the conductance oscillations, which contain additional information on the interaction of 
the charge carriers with the defect, requires the solution of the scattering problem of the 
influence of subsurface defects on the conductance of a small tunnel contact. 

In this paper we review a series of publications in which the theory of the electronic 
transport though a tunnel point-contact in the presence of a single defect below metal 
surface was developed. The organization of this paper is as follows. The model of the tunnel 
contact and the basic equations that describe the effect of subsurface defects on the STM 
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conductance are presented in Sec. II. The solution of the Schrodinger equation for elections 
that tunnel through the contact and are scattered by the defect is given. In Sec. Ill a 
method to determine the defect positions below a metal surface is formulated on the basis 
of an investigation of the nonlinear conductance of the contact. A signature of the Fermi 
surface anisotropy in STM conductance in the presence of subsurface defects is discussed 
in Sec. IV. In Sec. V we present the results of investigations of the effect of a subsurface 
magnetic defect on the tunnel current, including the signature of a Kondo impurity and that 
of a magnetic cluster having an unscreened magnetic moment. In Sec. VI it is shown that a 
strong magnetic field leads to specific magneto-quantum oscillation periods which depend on 
the distance between the contact and the defect. The possibilities of studying the interference 
of quasiparticles in a superconductor is analyzed in Sec. VII. In Sec. VIII we conclude by 
discussing the possibilities for exploiting these theoretical results for sub-surface imaging 
along with experimental investigations of physical characteristics of subsurface defects. 



II. QUANTUM INTERFERENCE OF SCATTERED ELECTRON WAVES IN 
THE VICINITY OF A POINT CONTACT 



A. Model of a STM contact, and the Schrodinger equation for the system 



As a model for the STM experiments we choose an inhomogeneous tunnel contact between 
two metal half-spaces separated by an infinitely thin interface. The potential barrier in the 
plane of the interface, at z = 0, is taken to be described by a delta function [36], 

U(t) = U f(p)6(z), (1) 

where p is the radius vector in the plane of the interface, perpendicular to the z axis. The 
function f(p) — > oo at all points of the plane z = except for a small region defining the 
contact, having a characteristic radius a, at which f{p) is of order 1. As an example, a 
suitable model for the function f(p) for the "STM tip" is the Gaussian function f(p) = 
exp(p 2 /a 2 ) with small a. Another useful model of the junction is an orifice of radius a for 
which f(p) = 1 for p < a in the plane of the contact (Fig. [[]). 

Of course, such a model describes only the qualitative features of the conductance of 
an STM contact, and does not contain such parameters as the tip radius, or the distance 
between the STM tip and the sample as is represented, for example, in the model by Tersoff 



and Hamann |43] . In principle these properties of the system may be included in the model 
as parameters of the function f(p). The advantage of the model by Kulik et al. [36[ is 
the possibility of finding exact analytical solutions of the Schrodinger equation in the limit 
Uq — > oo. The equations are considerably simplified in the case of a small contact a — > 0. 
The wave functions obtained in the framework of the model barrier ([1]) properly describe 
the spreading of electron waves into the bulk metal from a small region on its surface. A 
numerical value for the STM conductance plays the role of a scale factor for the conductance 
oscillations, and for the further considerations below it is of less importance. 
A defect in the vicinity of the interface can be described by the potential 

D(r) = gD (\r-r \), (2) 

where g is the constant of the electron interaction with the defect, and Do(\r — tq\) is a 
spherically symmetric function localized within a region of characteristic radius ro centered 
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at the point r = ro, which satisfies the normalization condition 

4tt I dr'r' 2 D (r') = 1. (3) 



The electron wave function ip (r) in a metal with a dispersion relation e (k) must be found 
from the Schrodinger equation [44| 

e(k- J: A) + ag e fx B H + eV{r) +D{r) + U{r) $ = eifj. (4) 

Here k = — iV, A(r) is the vector-potential of the stationary magnetic field H, and V(r) 
is the applied electrical potential, a = ±1 corresponds to different spin directions, /i B = 
ehjlrriQC is the Bohr magneton, where m is the free electron mass, and g e is the electron g- 
factor. The function ip (r) satisfies at z = the following boundary condition for continuity 
of the wave function 

V>(p,+0)=V(p,-0), (5) 

and the condition, which for a 5— function barrier is obtained by the integration of the 
Schrodinger equation fl3j) over an infinitesimal interval near the point z = 

+o 

J dzs (k - J-a) V (P, ^) = -Uof(p)^ (p, 0) . (6) 

-0 

In this section below we consider a solution of Schrodinger equation (j3j) for a free electron 
model with an electron effective mass m* and a dispersion relation e (k) = h 2 k. 2 /2m* in the 
absence of external fields (H — 0,V — ). In this case the condition (j6]) reduces to the 
well-known condition for the jump of the derivative of the wave function 

ft (P, +0) - V4 ("P, 0) = -^f (p) V (P 0) • (7) 

The effects of applied voltage, Fermi surface anisotropy and magnetic field are discussed in 
next sections. 



B. Wave function due to an inhomogeneous tunnel barrier 

Here we follow the procedure for the finding the electron wave function in the limit 



Uq — > oo that was proposed in Ref. [36|. To first approximation in the small parameter 1/Uq 



the wave function ip (r) can be written as: 

i) (r) = % (r) + y? (r) , (8) 

where <p is of order 1/Uq. This latter part of the wave function (JSJ) describes the electron 
tunnelling through the barrier and determines the electrical current. The first term in the 
Eq.(JHD is the solution of the Schrodinger equation for the metallic half-spaces without the 
contact 

^ ( r ) = e iKp [e kM - e~ tkzlzl ) , (9) 
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where n and k z are the components of the wave vector k parallel and perpendicular to the 
interface, respectively. The expression fl9]) satisfies the boundary condition ip Q (p, 0) = at 
the interface. 

Substituting the wave function (|S]) into the boundary conditions ©and (jBD one must 
match terms of the same order in 1/U . As a result the conditions (EJ), (EJ) are reduced to 

<p o (p,+0) = <p o (P,-0), (10) 

t (k z ) e iKp = f (p) <p (p, 0) . (11) 

where 

t (k z ) = h 2 k z /im*U ] \t\ < 1, (12) 

is the amplitude of the electron wave function passing through the homogeneous barrier. 
Developing the function <p (p, z) as a Fourier integral in the coordinate p, and using the 
Eq. dTTJ, we find H 

oo oo ; 

^ (p, z ^ 0) = ^ | d^'^N | ^'^7^, (13) 



-oo 



where fc^ = \/k 2 — k' 2 . For a homogeneous 5-function barrier, /(p) = 1, Eq. (fl3|) transforms 
into a transmitted plane wave having an amplitude t. 

The characteristic radius of the region on the surface through which electrons tunnel 
from the STM tip into the sample is of atomic size, a ~ O.lA, while the Fermi wave vector 
k F ~ lA" 1 . By using the condition fcp a i,2 1 after integrating over k,' and p' in Eq. fll3j) 
we find 

^oW^^^f'W. (14) 

The incident plane wave is transformed into a spherical p-wave hi (kr) tHM after scattering 
by the point contact. In Eq. (JHJ), and below, (x) are the spherical Hankel functions. 
Note that the wave function <Po( r ) P^P ^ s zero i n au poi^s on the surface z — 0, except the 
point r = (at the contact) where it diverges. This divergence is the result of taking the 
limit a — > in the integral expressions for <Po( r ) ([I3D - Yet, Eq. (|14p gives a finite value for 
the total charge current through the contact as obtained by integration over a half-sphere 
of radius r with its center in the point r = for r — > 0. 



C. Electron scattering by a single defect in the vicinity of a tunnel point contact 

As a result of current spreading only a small region near the point contact noticeably 
influences the conductance. For high purity samples only a few defects will be found in 
this region. At low temperatures the distance between the contact and the nearest defect, 
ro, is smaller than the electron mean free path due to electron-phonon scattering and the 
electrons are elastically scattered by the single defect only. The wave function of transmitted 
electrons, <p (r), which takes into account the scattering by the defect, can be expressed in 
terms of the retarded Green function G$ (r, r'; e) of the homogeneous equation ([!]) at D — 0, 
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FIG. 1: Model of the tunnel point contact as an orifice in an interface that is nontransparent for 
electrons except for a circular hole, where tunnelling is allowed. Trajectories are shown schemat- 
ically for electrons that are reflected from or transmitted through the contact and then scattered 
back by a defect. 



in absence of impurity scattering. To first approximation in the transmission amplitude t 
( jl"2l) the integral equation for (p (r) is given by 

<p(r) =<p (r) +g J dr'D (|r' - r |) G+ (r, r'; e) <p (r') , (15) 



{hP(k\T-T>\)-hW(k\T-T'\)}, (16) 



where 

is the electron Green's function of Eq. (j4j) for the semi- infinite half-space (Uq — > oo), 
r' = (p', — z') ,and y? ( r ) is given by Eq. f lT3|) . For small g Eq. (fl5j) can be solved by 
perturbation theory, i.e. in first approximation in g the function if (r') in the integral term 
should be replaced by (p (r') ■ 

For a short range potential (Atd <C 1) the function <£> (r') can be taken outside of integral 
in Eq. (|T5|) and the scattered wave function is written as [45] 

<p{r) = <p {r) + T (k) <p (r ) G+ (r, r ; e) , (17) 

where ^ 

T W = l-p/dr'Do(|r'-r |)G+(r ,r';e)" (18) 

Note that Eq.lfPTj) is valid far from the defect (|r — r | 3> ro) and the function D (|r' — r |) 
must provide the convergence for the integral in the denominator of Eq. (fT8|) at r' — > tq. As 
is well known, s-wave scattering is dominant for scattering by a short range potential and 
the scattering matrix (TT8]) can be expressed by the s-wave phase shift 5 [46j 



ink 2 e 2iSo 
™*k i + i ( e 2i5 _ i) fr« ( 2 foo)' 



r (*) = E=r : rr^ zz^r, ■ (19) 
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FIG. 2: Spatial distribution of the square modulus of the wave function in the vicinity of the 
contact in the plane perpendicular to the interface pas sing through the contact and the defect. 
Distances are given in units of the inverse wave number 51]. 



The effective T-matrix is an oscillatory function of the distance z$ between the defect and 
the interface that results from repeated electron scattering by the defect after its reflections 
from the interface. Figure [2] illustrates the spatial distribution of the square modulus of the 
wave function ( TT71) in the vicinity of the contact with the defect placed at ro = (5, 0, 15) /k. 

III. FRIEDEL-LIKE OSCILLATIONS OF THE TUNNEL POINT CONTACT 
CONDUCTANCE 

A. Voltage dependence of the STM conductance 



In the case of a small transparency ( 1121) the applied voltage drops entirely over the 
barrier and the electrical potential can be chosen as a step function V^(r) = VQ(—z). At 
zero temperature electrons tunnel to the lower half-space when eV > 0, and for eV < 
electrons can tunnel only to available states in the upper half-space (Fig. [1]). 

The tunnelling current /(V) = I^ + \V) — I^~\V) is the difference between two currents 
flowing through the contact in opposite directions. Each of them can be evaluated by means 
of the probability current density (V) integrated over a plane z = const, and integrating 
over all directions of the electron wave vector 



4 ±) (V) = u(e) I dpQ{±z)(e(±k z )Re<p*(r)vMr)) 



e ' 



(20) 



where v (e) is the electron density of states for one spin direction, (...) E denotes the average 
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over an iso-energy surface e (k) = e, 



(..->« 



/ 



V(k)= 



dSu 



(21) 



e(k)=e 



dSk is an element of the .so-energy 8Ur faee in k-space, and v = ^ is the veloc.ty operator. 
For a free-electron model of the energy spectrum v z = -^r Jj- A voltage dependence of the 
current density J^(V) ( 1201) is defined by the dependence of the absolute value of the wave 
vector for the incident on the contact electron |k (V) \ = a/ k 2 — 2m* \eV\ /h 2 . 
The total current through the contact is 



I(V) 



<7=±1 

r(-) 



ji + \V)f F (e-eV)(l-f F (e))- 



Ji'(V)f F (e)(l-f F (e-eV)) 



(22) 



where /p (e) is the Fermi function. 

The current- volt age characteristics / (V) is calculated by substituting wave function ( TT71) 
into Eq. (122]) and taking into account Eqs. (|T3|) and (|T6|) . Retaining only terms to first order 
in g (i.e. ignoring multiple scattering at the impurity site, in Eq. ( ITT)) T(k) ~ g), and in the 



limit of low temperatures, T = , the conductance G(V) = dl/dV can be written as [51 



G(v ,V) 



47r 3 (m*U ) 



dp x dp 2 



/(Pi)/(p 2 ) 



x 



(23) 



fc F fc F^fc F (Pl,P 2 ) 



2 / k 5 dkF k { Pl ,p 2 ) 



where 



^fc (pi,P 2 ) 



h (kp) 
kp 



Am*gkj l (kp) z\ 
nh 2 kp A1A2 



ji (k\i) yx (k\i 



(24) 



|Pi - P2I , 



iPo -P 
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A, 



Zq + \p — p 2 \ 2 , and ji(x) and yi(x) are the 



spherical Bessel functions, and 



k F (V) = Jk 2 + 2m*eV/h 2 



(25) 



is the Fermi wave vector k F accelerated by the potential difference. In Eq.( l23|) for definiteness 
a positive sign of the bias is chosen, eV > 0. 

If the contact radius a < Ap (Ap = l/k F is the Fermi wave length), the expression for 
the conductance (1251) can be simplified 



rl W k F r 



v (k F r ) 



(26) 
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where 



Go 



I* (Ml 



2 e 2 (k F a) 4 
36nh 



is the inherent conductance of the tunnel point contact, r$ 

1 



iPol 



q(x) 



l + x X 



w (x) = — [(x 2 — l) sin 2x + 2x cos 2x\ 



v (X) 



X ' 



[2x (Ax 2 - 7) sin 2x + (2x A - Ux 2 + 7) cos2x] 



and 
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6m* fo 



71 



h 2 



-9, 



(27) 



(28) 

(29) 
(30) 

(31) 



is the dimensionless constant of interaction. 

Equation ( 126]) describes the oscillations of the STM conductance as a function of the 
distance r between the STM tip and the subsurface defect, and as a function of the bias 
eV. For distances between the contact and the defect r 3> Ap and eV <C the oscillatory 
dependence becomes sinusoidal 



G(r ,VO-G cx-£sin2W 



(32) 



Oscillations of the STM conductance as a function of the voltage, due to the quantum 
interference caused by impurity scattering, were observed by Untiedt et al. [48] , and Ludoph 
et al. Hi . 



B. Determination of the defect positions 

Now we proceed to discuss whether this effect can be exploited experimentally for three 
dimensional mapping of subsurface impurities. The position of the defect in the plane parallel 
to the surface can be found from an analysis of oscillatory pattern in the dependence G (p ) . 
In the majority of cases the center of this pattern corresponds to the tip position directly 
above the defect, p = 0. A possible effect of the Fermi surface anisotropy is discussed in 
the next section. Note that, in contrast to the case of surface defects, the oscillations in the 
conductance (1261) are not periodic in the tip distance p along the surface, but their period 
is defined by the distance r = a/Po + z o- Generally, the depth z may be found by fitting 
the experimental data to the theoretical dependence G (p , z ) (T26|) . Figure 3 illustrates the 
oscillatory component of G (r , V — 0) as a function of p for different choices of z . In this 
plot we have used the values for the constant g — 1, the Fermi wave vector k F = 1.360 
A -1 and the interatomic distance d = 1.805 A for Cu. Thus, the plots correspond to defect 
positions in the third, fourth, and fifth layers below the Cu surface. The dependencies G (p) 



closely resemble the observations by Quaas et al. [32| for Co atoms embedded in Cu(lll). 

For the determination of the defect depth zq one may use the periodicity in phase d = 
2kpr of G(r , V) (1321) at sufficiently large r . According to Eq. (1321) at V = two sequential 
radii p 01 and p 02 , p 02 > p i> corresponding to neighboring maxima (or minima) satisfy the 
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obvious condition of periodicity, Ai? = 2k F (^/ p\ 2 + z\ — a/PoT+^o) = For known k F 
it is a simple algebraic equation for z®, the solution of which is 



z ° = TIT V fc F (P02 - Pgr) - 2vr 2 fc| (pg 2 + pg ) + ^ (33) 

Note that k F (p 02 — p 01 ) > n and the radicand is positive. A second possibility of changing 
the product k F ro is by varying the maximum value of the electron wave vector by the applied 
voltage. 

A first approach for determining the defect depth z from the bias dependence of the 
period of the Friedel-like oscillations of the STM conductance was described by Kobayashi 
[42|. The depth z can estimated by tracing the point p while changing the bias voltage \eV\, 
keeping the phase of the oscillations d constant: k F a/pq + z$ = k F (V) \J p'§ + Zq, where p 
and Pq are the positions corresponding to two different bias voltages V — > and V2 = V 
but the same phase (for example, a fixed maximum) [42]. The solution of mentioned above 
equation with kp (V) (J2"5j) gives z 



Zo = ^M-g-^? , (34) 

where eV > 0, p > p' . 

The method proposed in Ref. 51] has certain advantages. If the STM tip is placed above 
the defect (|p | z ) the conductance amplitude decreases with depth of the defect as z%, 
which gives hope to observe the defects at sufficiently large distances below the surface. 
The depth of an impurity may be derived from the G(V) curve at p = 0, which shows 
oscillations in eV with period eAV, and 

z = ^ . 35 

k F - k F (Ay) 

In a real experiment it is not necessary to observe a full period of G(V) and, for example, 
a quarter of the period will be sufficient for the determination of the defect depth [5l[ . 



IV. SIGNATURE OF THE FERMI SURFACE ANISOTROPY 

In most metals the dispersion relation for the charge carriers is a complicated anisotropic 



function of momentum. This leads to anisotropy of the various kinetic characteristics (44 
Particularly, as shown in Ref. [50[ , the current spreading may be strongly anisotropic in the 
vicinity of a point-contact. This effect influences the way the point-contact conductance 
depends on the position of the defect. For example, in the case of a Au(lll) surface the 
'necks' in the Fermi surface (FS) should cause a defect to be invisible when probed exactly 
from above. 

Qualitatively, the wave function of electrons injected by a point contact for arbitrary FS 
e (k) = e F has been analyzed by A. Kosevich [50|. He noted that at large distances from 
the contact the electron wave function for a certain direction r is defined by those points on 
the FS for which the electron group velocity is parallel to r. Unless the entire FS is convex 
there are several such points. The amplitude of the wave function depends on the Gaussian 
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FIG. 3: Dependence of the normalized oscillatory part of the conductance on the STM tip position 
for different depths zq of the defect below the surface; ro = (p , 0, z a ) ,g = 1. 

curvature K in these points, which can be convex (K > 0) or concave (K < 0). The parts of 
the FS having different signs of curvature are separated by lines of K = (inflection lines). 
In general there is a continuous set of electron wave vectors for which K = 0. The electron 
flux in the directions having zero Gaussian curvature exceeds the flux in other directions 

Electron scattering by defects in metals with an arbitrary FS can be strongly anisotropic 
. Generally, the wave function of the electrons scattered by the defect consists of several 
superimposed waves, which travel with different velocities. In the case of an open FS there 
are directions along which the electrons cannot move at all. Scattering events along those 
directions occur only if the electron is transferred to a different sheet of the FS jiij . 

In this section we analyze the effect of anisotropy of the FS to the signals for determination 
of the position of a defect below a metal surface by use of a STM. We show below that the 
amplitude and the period of the conductance oscillations are defined by the local geometry 
of the FS, namely by those points for which the electron group velocity is directed along the 
radius vector from the contact to the defect. General results are illustrated for the FS of 
noble metals. 

At first we do not specify the specific form of the dependence e (k), except that it satisfies 
the general condition of point symmetry e(k) = e (— k). In the reduced zone scheme a 
given vector k identifies a single point within the first Brillouin zone. As for isotropic FS, 
the electron wave function ip (r) in the metal with an arbitrary dispersion relation can be 
found at U — > oo by using the method described in Sec. III. The boundary conditions for 
the transmitted wave function ip (r) have the same form as Eqs. fllOp . fill I) in which the 
function t (k) must be replaced by 

(36) 

z=-0 

For the model of free electrons Eq. fl36|) transforms into Eq. (]12p . The components of vector 
k perpendicular to the interface for electrons incident on the contact, k z (k, e), and reflected 
from the contact, k T z ei (k,e), are related by conditions of conservation of the energy e and 
the tangential component n of the wave vector 

e fV n , k z ) = e (V ef , kf) = e; rc in = K rcf = k. (37) 



ioz 



,ik x \x\ 



-ik r z et \z\ 
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The wave function scattered by the defect is defined by the general relation ( 1X5]) . 

General expressions for the STM conductance into a metal having an arbitrary FS one 
can find in Ref. [55] . Here, we present simplified asymptotic expressions for the oscillatory 
part of the conductance AG^* (r ) (the difference between the total conductance and its 
value in the absence of the defect) which are valid for large distances between the contact 
and the defect , r ^> Ap [55l . 



E 



AGS (ro) = M (I* (k)| 2 e K)> ep ■ (38) 

1 

: sin(/i (k 0s ) r + <j> 8 ) cos (h (k 0s /) r + <f) s ,) . 



VI^(ko a )^(k 



Os' 



All functions of the wave vector in Eq. ( 138]) are taken at the points of the FS for which the 
electron group velocity v is parallel to the vector r = r n , h (e F , k ) = k n , k is the 
wave vector corresponding to the point on the FS, in which Vo||n . The function /i(£p,ko) 
is well known in the differential geometry as the support function of the surface e (k) = e-p 
13]. If the curvature of the FS changes sign, there is more than one point k 0s (s = 1,2...) 



for which vo s || no- It may also occur that for given directions of the vector ro v no for all 
points on the FS, and the electrons cannot propagate along these directions [441 ] . For such 
ro the oscillatory part of the conductance is zero. 
In Eq. dSHD 

T^M^T (l + sgntf(ko)), (39) 



4 ° \ dk 

..(+) _ l(+) 



X 



(. . . ) ep is defined by Eq. 021]), W ' = kY 1 (k x , k y , £p) in the point defined by the direction 
of the vector n in k-space, K (k ) ^ is the Gaussian curvature of the FS, 

K (e F , k ) = - — ^ y~] A ik n 0i noj , (40) 

Vo • 

i u i i,j=x,y,z 

Aij = 8 ^-™ k ) is the algebraic adjunct of the element m^ 1 (k) = p of the inverse 
mass matrix m _1 . 

The Eq. (|38p is valid, if curvature 7^ 0. For those points at which K the amplitude of 
the electron wave function in a direction of zero Gaussian curvature is larger than for other 
directions. This results in an enhanced current flow near the cone surface defined by the 



condition K = [50|, |55|. If the FS is open, there are directions along which the electron 
flow is absent. These properties of the wave function manifest itself in an oscillatory part of 
the conductance (138]) : 1) The amplitude of oscillations is maximal if the direction from the 
contact to the defect corresponds to the electron velocity belonging to an inflection line. 2) 
There are no conductance oscillations, AG*J = 0, if this direction belong to cones, in which 
the electron motion is forbidden. 

For an ellipsoidal FS the Schrodinger equation can, in fact, be solved exactly in the limit 
a — > 0, Uq — >■ 00 and the conductance of the contact can be found for arbitrary distances 
between the contact and the defect. For this FS the dependence of the electron energy e on 
the wave vector k is given by relation, 

«M = t E g; («) 
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were ki are the components of the electron wave vector k, 1 / are constants representing 
the components of the inverse effective mass tensor m _1 . 

Accurate to within first order in q (i.e. ignoring multiple scattering at the impurity site), 



the conductance in the limit V — > 55] is given by 

6g(2e F f 2 ( z \ 2 



G cl1 (r ) = Gf 



w (h (ko) r ) 



7rh 5 y/m zz det [m x ] \/i(ko)r 0/ 
where Gq 11 in Eq. (142[) is the conductance in the absence of a defect (g = 0) |55j: 



(42) 



G cii = ^ae F . (43) 

9irh 3 Ugy/m Z2 det [m- 1 ] 



1/2 

(44) 



w (At) is given by Eq. fj2U|) . 

The center of the oscillation pattern in the conductance G el1 (ro) as the function of the 
tip position p corresponds to p = p 00 with respect to the point contact at r = 0, where 

/ m zz m zz \ 

Poo = z o , • 45 

\m zx m zy J 

The support function h for such tip position, 



k n 00 = k zF = -\/2e F m zz ; (46) 

corresponds to the extremal value of the chord 2k zF of the FS in the direction normal to the 
interface, n 00 is the unit vector in the direction of the vector r o = (Pocb ^o)- Figure H] shows 
that AGq] c = G cl1 — Gq 11 is an oscillatory function of the defect position p that reflects the 
ellipsoidal form of the FS and the oscillations are largest when the contact is placed in the 
position p 00 , defined by Eq. f T45|) . 

In deriving Eq. (l38|) it has been assumed that eV — > 0. For finite voltage, but eV e F , 
all functions of the energy e in Eq. ( 1381) can be taken at e = e F , except h (e, k ) in the 
oscillatory functions. When eV e F , 

Oh Oh cV 
h{e F + eV,k )^h(e F ,k ) + —eV, k F — (47) 

OB F u£ F B F 

and when the product (eV/e F )k F ro 3> 1 clearly the conductance (1381) is an oscillatory func- 
tion of the voltage V. The periods of the oscillations are defined by the energy dependence 
of the function h(e,k ). The results obtained properly describe the total conductance at 
eV e F and also can be used for the analysis of the periods of the oscillations at eV < e F . 

Further calculations require information about the actual shape of the FS, e (k) = e F . In 
Ref. [55j a model FS in the form of a corrugated cylinder was considered. Using this model, 
for which analytical dependencies of the conductance on defect position can be found, the 
manifestation of common features of FS geometries in the conductance oscillations was 
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FIG. 4: Dependence of the oscillatory part of the conductance, AG^ C , as a function of the position 
of the defect p i n the plane z = zq. The shape of the FS (liTj) is defined by the mass ratios 
m x /Tn z =l, m y /m z = 3, and the long axis of the ellipsoid is rotated by 7r/4 around the x— axis, 
away from the y— axis. The coordinates are measured in units l/k z p (|46p and the defect sits at 
z = 5 @. 

described: the anisotropy of convex parts ('bellies'), the changing in sign of the curvature 
(inflection lines), and the presence of open directions ('necks'). 

In Ref. (68[ a numerical analysis of the conductance oscillation pattern was made for the 
noble metals copper, silver and gold on the basis of Eq. ( 1381) . The parameterization of the 
FS was taken from 1561, 



e(k) = a 



—3 + cos cos— — h cos— ^ cos V cos cos h (48) 

2 2 2 2 2 2 y 1 

'—3 + cosk x a + coskyd + cosk z a)] , 



which is accurate up to 99%. The values for the constants are r = 0.0995, e/ot — 3.63, and 
a is different for each metal. For copper, silver, and gold a = 0.361nm, a = 0.408nm, and 
a = 0.407nm, respectively. The Fermi energy of copper is 7.00eV, for silver 5.49eV and for 
gold it is 5.53eV. 

The results of computations for three crystallographic orientations are presented in Fig. |5J 
All distances in Fig. [5] are given in units of Ap, which for copper is 0.46nm, and for silver and 
gold it is 0.52nm. For each of the surface orientations the graphs have the symmetries of 
that particular orientation of the FS. In all figures 'dead' regions can be seen, for which the 
conductance of the contact is equal to its value without the defect, showing no conductance 
oscillations. These regions originate from the 'necks' of the FS and their edges are defined 
by the inflection lines. For all orientations of the metal surface the defect position in the 
plane of the surface corresponds to a center of symmetry. The appearance of 'dead' regions 
depends on the depth of the defect, which can be estimated in the following way: The 
orientations of the 'neck' axes define the axes of the cones with an opening angle 27, in 
which there are no scattered electrons. Vertexes of the cones coincide with the defect. The 
radius R of the central 'dead' region, R = z$ tan (7) , is proportional to the depth of the 
defect (32l ] . 

The possibility of visualizing the Fermi surface of Cu in real space by investigation of 
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FIG. 5: a) Fermi surface described by Eq. (|48p relative to the contact axis for three principal lattice 
orientations, b) Gray scale plot of the tunnelling point contact conductance G as a function of the 
contact position for a defect at the origin, at a depth of 5Ap and for a (100) surface plane; the x and 
y directions each correspond to (100) directions, c) Same plot for a (111) surface orientation; the 
x and y directions correspond to [112]and [110] directions, respectively, d) Same plot for a (110) 
surface orientation; the x and y directions correspond to [001]and [110] directions, respectively 68]. 



interference patterns caused by subsurface Co atoms has been demonstrated by Weismann 
et al. 33 . 



V. SUBSURFACE MAGNETIC DEFECTS 



A. Kondo impurity 



In the case of a magnetic defect at low temperatures (T <C Tk, where Tk is the 
Kondo temperature) the Kondo resonance results in a dramatic enhancement of the ef- 
fective electron-impurity interaction [53| and perturbation methods become inapplicable. 
Kondo correlations give rise to a sharp resonance in the density of states at the energy 
e(k) = £k near the Fermi level. For e(k) — > £k the effective electron scattering cross section 



acquires a maximum value corresponding to the Kondo phase shift £ok = tt/2 53]. In this 
case, multiple scattering needs to be taken into account, even for a single defect, because of 
electron reflection by the metal surface. 

In this subsection the conductance is expressed by means of a s-wave scattering phase 
shift 5q. The results describe the influence to the conductance of multiple scattering of 
the electrons, which results in the appearance of harmonics in the dependencies of G on 
the applied voltage, and on the distance between the contact and the defect. The analysis 
of the non-monotonic voltage dependence of the conductance is applied s pec ifically to the 
interesting problem of Kondo scattering, using an appropriate phase shift [54|: 



8 (k) 



IT 

--tan 



-i 



e(k) - e K 



K 



+ 5, 



0D- 



(49) 



17 



The first term in Eq. (|49j) describes the resonant scattering on a Kondo impurity level Ek {Tk 
is the Kondo temperature). For a non-magnetic impurity this term is absent. The second 
term Sqd takes into account the usual potential scattering. 

Taking Eq. ( II 7p for the wave function of a spherical Fermi surface and Eq. (1191) for the 
scattering matrix makes it possible to find the differential conductance, G = dl/dV, of 
the tunnel point contact in the approximation of s-wave scattering. For \eV\ < e F and for 
eV > 0, G(V) is given by Q 



G(V) = G c 



q(V) + 




k 6 

n I 



dkk 5 <S> (k) 



kp 



(50) 



and for eV < 0, 



G(V) = G CJ 




4 J 



k 2 - 



2m*eV 

h 2 



$(fc) 



(51) 



Here Go is given by Eq. (|27|) . k F = \/k\ + 2m*eV/h 2 , and 



$ (k) = F- 1 sin 5 % [12J1 (Ar ) (-yi (Ax ) cos 5 + 



(52) 



{ji {kr ) (j (2kz ) 



y (2kz ) y 1 (kr )} sin 5 ) + 



6 (1 - jo (2fe8o)) (^o)" 4 (1 + (fcr ) 2 ) sin5 ] 



2 sin 5 x 
1 

2 (2fcz ) ; 



- jo(2/c2o) sin5 - yo(2kz ) cos5 c 



(53) 



#0 (fc) is s-wave phase shift SSI and ji(x) and yj(a;) are the spherical Bessel functions. 

At low voltage the conductance can be expressed by an expansion in the small parameter 
1/ (k F z ) < 1, 



f 2 i 00 

G(0) = G 1 + 12^|— -^(-1) 



sin n 5r, 



n-1 



(k F r ) 



r o (^f^o) ^ (2k F z 
2 t sin (2k F (r + (ra — 1) z ) + n5 



(54) 



/c F r 



cos (2fc F (r + (n — 1) 2 ) + n5 c 



The second term in Eq. ( 1541) gives the sum over n scattering events by the defect and n — 1 
reflections by the surface. If we keep only the term for n — 1 Eq. (I54j) reduces to the result 
obtained by perturbation theory in Sec. 3 above, which is valid for 5n ~ —gm*k F /2irh 2 <C 1. 
The arguments of the sine and cosine functions in Eq. (|54"|) correspond to the phase that the 
electron accumulates while moving along semiclassical trajectories. 
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FIG. 6: Difference 5Gk{V) / Go between the voltage dependencies of the conductance for a magnetic 
and a non-magnetic impurity. The parameters ek = 0.9ef, Tk = 0.01£f, and <5od = 0.1 in Eqj4"9l 
are used 1461. 



The voltage dependence of the conductance is not symmetric around V — 0. This asym- 
metry arises from the dependencies of the phase shift So{kp) (l4"9~j) and the absolute value of 
the wave vector k F = \Jk^. + 2m*eV/h 2 on the sign of eV. The physical origin of this asym- 
metry comes from the fact that the scattering amplitude depends on the electron energy in 
the lower half-space (see Fig. [T]), where the defect is situated. This energy is different for 
different directions of the current. 

It is interesting to observe that the sign of the Kondo anomaly depends on the distance 
between the contact and the defect Tq. This distance in combination with the value of the 
wave vector kp determines the period of oscillation of G(V). If the bias eV^ coincides with 
a maximum in the oscillatory part of conductance the sign of the Kondo anomaly is positive 
and vice versa, a negative sign of the Kondo anomaly is found at a minimum in the periodic 
variation of G(V). 

The FigJB] shows the difference SGk(V)/G = {G m — G n )/G between voltage dependen- 
cies for a magnetic G m and a non-magnetic G n impurity, having the same potential scattering 
strength. The plots in Figj6] demonstrate the evolution of the shape of the Kondo anomaly 
for several values of the distance between the contact and the impurity, placed on the contact 
axis. The change of distance changes the periodicity of the normal-scattering oscillations 
which leads to a changing of sign in the Kondo signal. A similar dependence of the differ- 
ential conductance with the distance between an STM tip and an adatom on the surface of 
a metal has been obtained theoretically in Refs. 57, 58[ in terms of a Anderson impurity 
Hamiltonian 5^]. Note that we obtain a Fano-like shape of the Kondo resonance in the 



framework of a single-electron approximation |46|, while in Refs. [57|, [58| many-body effects 
were taken into account. 



B. Magnetic cluster 

In this subsection we consider the influence on the conductance of a tunnel point-contact 
between magnetic and non-magnetic metals of a defect having an unscreened magnetic 
moment, in a spin-polarized scanning tunnelling microscope (SP-STM) geometry [24]]. A 
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magnetic cluster is assumed to be embedded in a non-magnetic metal in the vicinity of 



the contact. As first predicted by Frenkel and Dorfman |61[] particles of a ferromagnetic 
material are expected to organize into a single magnetic domain below a critical particle size 
(a typical value for this critical size for Co is about 35nm). Depending on the size and the 
material, the magnetic moments of such particles can be fi cS ~ 10 2 — 10 5 /i B [62]. 

Generally, the moment /Lt eff of the cluster in a non-magnetic metal is free to choose an 
arbitrary direction. This direction can be held fixed by an external magnetic field H, the 
value of which is estimated as H ~ T/fi cS , where T is the temperature (see, for example, 



Ref. [62|). For /i eff ~ 10 2 /i B and T ~ IK the field H is of the order of 0.01T. If H is much 
larger than the magnetocrystalline anisotropy field of the magnetic STM tip, the direction 
of the external magnetic field controls the direction of the cluster magnetic moment but 
its influence on the spin-polarization of the tunnel current is negligible. In this case the 
magnetic moment fj, eS of the cluster is 'frozen' by the field H and the problem becomes a 
stationary one. 

If the external magnetic field is sufficiently weak and the radius of the electron trajectories 
7h = hck-p/eH is much larger than the distance between the contact and the cluster r , the 



effects of modulation of the tunnel current due to electron spin precession [69| and trajectory 



magnetic effects [66j are negligible. 

The geometry of a SP-STM experiment can be described in the framework of the model 
presented in Fig. [TJ in which the half-space z < is taken up by a ferromagnetic conductor 



with magnetization M. In Ref. [67[ the direction of the vector M , which defines the 
direction of the polarization of tunnel current is chosen along the z-axis. In real SP-STM 
the polarization of STM current is defined by the magnetization of the last atom of the tip 
[24j . A magnetization oriented along the contact axis can be obtained, for example, for a 
Fe/Gd-coated W STM tip 0. 

The interaction potential D(r) of the electrons with the cluster is a matrix consisting of 
two parts 

1 



D(r) = Igl + —J » cS *) D (\t - r |), (55) 

where g is the constant describing the non-magnetic part of the interaction (for g > the 
potential is repulsive), J is the constant of exchange interaction, /x off = yU eff (sina, 0, cos a) 
is the magnetic moment of the cluster, a = (a x ,a y ,a z ) with a M the Pauli matrixes, and / 
is the unit matrix. The function D(r) satisfies the condition Q. In the case of spin flip 
scattering the spinor electron wave functions satisfy the Schrodinger equation in which 
the scattering potential must be replaced by the matrix D(r) ( 153]) . Under the assumptions 
that the potential D(r) and the transparency of the tunnel barrier in the contact plane are 
small the two-component wave function can be found by the method described in Sec. II. 

The difference in absolute values of the wave vectors k CT for spin-up and spin-down elec- 
trons (for the same energy e), which move towards the contact from the ferromagnetic bank, 

k n = jrv^m* (e=F4ng ef i B M), (56) 

results in different amplitudes t a = t(k a ) (see Eq. ([12]) ) of the electron waves injected into the 
non- magnetic metal for different directions of the spin (g e is the electron g-factor). The total 
effective polarization P eS of the current depends on the difference between the probabilities 
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of tunnelling for different a, 



Pe S (e) = KJS - (57) 
Ft I + \H\ 



The conductance G of the contact at T = and eV <C Ep is given by [6j 

(58) 

- £=ep 

where Go is the conductance of the contact in absence of the cluster 



G — — — Gq 



1 + — — g + - — P cff (ep) J cos a W(r ) 
irh 2 V 2 ^b / 



g » = («f + *w "T^S ' (59 » 

k-pa is the absolute value of the Fermi wave vector in the magnetic metal for spin direction 
a (see Eq. (|56|) ). and 



W(r ) = J dr'D (\r'-r \) (£\ w(kr'). 



(60) 

The function w is defined by Eq. ( 1291) . When the radius of action of the function 
Dq(\t — ro I) is much smaller than the distance between the contact and the center of the 
cluster, ro, W(tq) is an oscillatory function of kr$ for krp > 1, as for point defect with 
kr D 1 (see, Eq. f f26|) at V = 0), but the oscillation amplitude is reduced as a result of 
superposition of waves scattered by different points of the cluster. The integral W(r ) fl60|) 
can be calculated asymptotically for r ^> r D , kr Q ^> 1, and kr D > 1. For a homogeneous 
spherical potential -D (l r l) = ^ / o 1 @( r D — r ) {Yd is the cluster volume) the function W(r ) 
takes the form 

W(r )^3( Z ^) 2 sin2k yff , (61) 
V ; V^o; (2A:r ) 2 kd K ' 

where d = 2r D is the cluster diameter. The last factor in Eq. f lBTj) describes the quantum 
size effect related with electron reflections by the cluster boundary. Such oscillations may 
exist, if the cluster boundary is sharp on the scale of the electron wave length. Fig. [7] shows 
the dependence of the amplitude of the conductance oscillations on the cluster diameter. It 
demonstrates that a 7r-phase shift may occur resulting from interference of electron waves 
over a distance of the cluster diameter. 

In Eq. (158]) the term proportional to P e s takes into account the difference in the prob- 
abilities of scattering of electrons with different a by the localized magnetic moment /-i cS . 
It depends on the angle a between the tip magnetization and /x eff , as cos a. The same 
dependence was first predicted for a tunnel junction between ferromagnets for which the 



magnetization vectors are misaligned by an angle a [631 ]. and this was observed in SP-STM 
experiments [241 ]. 

Note that once the spin-polarized current-induced torque pulls the magnetic moment 
away from alignment with H, the cluster moment will start precessing around the field axis. 
The Larmor frequency is defined by the magnetic field due to combining the external field 
H and the effective magnetic field produced by the polarized current. The precession of 
the cluster magnetic moment gives rise to a time modulation of the SP-STM current as for 
clusters on a sample surface 25, 2^. 
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FIG. 7: Dependence of the the oscillatory part of the conductance on the tip position on the 
metal surface for a subsurface magnetic cluster with different cluster diameters. The /9 -coordinate 
is measured from the point p = at which the contact is situated directly above the cluster; 
r = (0,0,10)/**; ?=0.5; J 



m* fcp 
Mb ft 2 



J/i eff = 2.5; P eff = 0.4; a = (671]. 



VI. MAGNETO-QUANTUM OSCILLATIONS 

A. Conductance oscillations in perpendicular magnetic field 

In a strong magnetic field the STM conductance exhibits characteristic oscillations in 
magnetic field, which are attributed to Landau quantization. This effect has been observed 



in Ref. [64J and the energy dependence of the effective electron mass was determined. An in- 
fluence of the magnetic field on the interference pattern, which is produced by two adatoms, 
in the STM conductance has been investigated theoretically (65| and horizontal stripes re- 
lated to the Aharonov-Bohm effect were predicted. 

In Sec. Ill it is demonstrated that the dependence G(r , V) undergoes oscillations in r 
and eV resulting from the variation of the phase shift between transmitted and scattered 
electron waves. Here we discuss another way to control the phase shift between the interfering 
waves: an applied external magnetic field H produces oscillations of the conductance as a 
function of H. 

Let us consider the contact described in Sec. II, now placed in a magnetic field directed 
along the contact axis, H = (0,0, if). Figure [8] shows schematically the trajectories of the 
electrons that are injected into the metal and interact with the defect. 

In what follows the Schrodinger equation is solved along the same lines as in Sec. II, and as 
zeroth approximation we use the well-known wave function for an electron in a homogeneous 
magnetic field. In Ref. 66[ the dependence of the STM conductance on magnetic field 



has been obtained under the assumptions that the contact diameter a is much smaller 
than the magnetic quantum length, an = yJh/m*VL, the radius of the electron trajectory, 
ru = hkp/m,*Q, is much smaller than the mean free path of the electrons, I 3> ro, and 
the separation between the magnetic quantum levels, the Landau levels, HQ is larger than 
the temperature k^T , (Q = eH/m*c is the Larmor frequency). Although these conditions 
restrict the possibilities for observing the oscillations severely, all conditions can be realized, 
e.g., in single crystals of semimetals (Bi, Sb and their ordered alloys) where the electron 
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FIG. 8: Schematic representation of the electron trajectories in a vicinity of a point contact in an 
external magnetic field oriented along the contact axis. 



mean free path can be up to millimeters and the Fermi wave length Af ~ 10 _8 m. Under 
condition of the inequalities listed the dependence of the conductance of the tunnel point 
contact on H is given by j66| 



G(H) = G C (H) 



gm 



2tt 3 (jVpt + N n ) Va% 



n=0 



n'=0 



Here 



Xa( n > r o) = exp ( j L„(f ) exp ( ^z ^/2m* (e F + afi B H - e n ] 



(62) 
(63) 



£o = Po/Zoh, and L n (£) are Laguerre polynomials, e n = HQ (n + |), a = ±1 is the spin 
index, N-p a is the number of electron states for one spin direction per unit volume at the 
Fermi energy, 

2 lei H nmax 

Nf ° = J k~2 ^ 2m * (^F + crfi B H-e n ), (64) 

n max = [|^] is the maximum value of the quantum number n for which e n < ep, and [x] is 
the integer part of the number x, G c is the conductance in absence of a defect, 



G C (H) = {irk) 



, , ea 2 (iV Ft + iV F; ) x 2 



m*Un 



(65) 



The conductance (1651) undergoes oscillations having the periodicity of the de Haas-van 
Alphen effect that originates from the step-wise dependence of the number of electron states 
Nfo- flM]) on the magnetic field. At n max (e F ) S> 1, fi B H/ep <C 1 (semiclassical approxima- 
tion), Eq.( 155]) can be expanded in the small parameter HQ/sp 



G C {H) ~ G 



9 /MA 3/2 ^ _i s , £p ^ 
2\£f/ ^ (2s) 3/2 V M] 4 



(66) 



where Go is the conductance of the contact at H = (see Eq. fl59|) ). 
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FIG. 9: Oscillatory part of the conductance of a tunneling point contact with a single defect placed 
at kpp = 50, kpZQ = 30. The full curve is a plot for Eq. ([62]) . while the dashed curve shows the 
component AG\ for the semiclassical approximation, Eq. (|68l) . The field scale is given in units 
l/k F r H ; g = 0.5. 
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FIG. 10: Dependence of the oscillatory part of the STM conductance on the tip position for different 
values of magnetic field, zq = 30/fep, g = 0.5. 

The oscillatory part of the conductance, AG(H) = G(H) — G C (H), which results from the 
electron scattering on the defect, is plotted in Fig. [9]for a defect placed at (p, z) = (50, 30) / kp. 
Figure [10] illustrates the dependence of the conductance (162]) on p coordinate of the defects 
for different H. The beating of the oscillation amplitude due to the difference of electron 
energies for different spin is seen at higher magnetic field. 

The dependence plotted in Fig. O G (H) , f]6"2"]) . contains oscillations with different periods. 
The semiclassical asymptotes at Ml <C e-p of the expression for the conductance Eq. (I6"2"]) 
allows us to explain the physical origin of these oscillations. By using the Poisson summation 
formula in Eq. (162]) the part of the conductance AG(H) related with the scattering by the 
defect can be written as a sum of two terms 

AG(H) = AG i + AG 2 , (67) 
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each of them describing conductance oscillations with different periods, which are discussed 
in more detail below. 

B. Effect of flux quantization through the trajectory of the scattered electrons 

The first term AGi(H, r ) in Eq. (IBTI) describes the long-period oscillations 

AG*! (H, r ) = -G g-^ sin ( 2k F r - 2ir^- ). (68) 

ft F r V *Q / 

where $o = 2irhc/e is the flux quantum. The flux, $ = HS pr , is given by the field lines 
penetrating the areas of the projections S pr on the plane z = of the trajectories of the 
electrons moving from the contact to the defect and back (see, the trajectory 2 in Fig. [8]). 
Such trajectories consist of two arcs, and there are a lot of trajectories with different S pr . 
As was shown in Refs. [66] among these trajectories the signal is dominated by the one that 
has minimal area given by S pv = 2S seg . Here S seg = r 2 (9 — sin 26) is the area of the segment 
formed by the chord of length p and the arc of radius r = sin 9, with 9 the angle between 
the vector r and the z-axis, sin 9 = p /r . Therefore, the oscillation AGi disappears when 
the defect sits on the contact axis, p = 0. Obviously, the origin of these oscillations lies in 
the curvature of the electron trajectories in a magnetic field. As we can see from Eq. ([6"8l the 
oscillations in the conductance AGi have a nature similar to the Aharonov-Bohm effect (the 
conductance undergoes oscillations with a period $/$o) an d are related to the quantization 
of the magnetic flux through the area enclosed by the electron trajectory. For illustration of 
this fact in Fig. [9] the full expression for the oscillatory part AG(H) of the conductance (the 
second term in Eq. (162!) ) is compared with the semiclassical approximation AGi(H, p , z ), 
Eq fU). 

For the observation of the Aharonov-Bohm-type oscillations the position p of the defect 
in the plane parallel to the interface must be smaller then r#, i.e. the defect must be situated 
inside the 'tube' of electron trajectories passing through the contact. At the same time the 
inequality p > an must hold in order that a magnetic flux quantum $o is enclosed by the 
area of the closed trajectory. 

C. Effect of longitudinal focusing of electrons on the defect by the magnetic field 

The short-period oscillations originate from the effect of the electron being focused by 
the magnetic field, and is described by the term AG^-ff, r ) in Eq. fIBTj) . At p = this 
term can be written as 

A G2W ^i Go «(|l) 3/2 £ <zl)! cos (, Fro + 2ls |L + _^). m 

In the absence of a magnetic field only those electrons that are scattered off the defect in the 
direction directly opposite to the incoming electrons can come back to the point-contact. 
When H ^ the electrons move along a spiral trajectory and may come back to the contact 
after scattering under a finite angle to the initial direction (the trajectory 1 in Fig.8). For 
example, if the defect is placed on the contact axis an electron moving from the contact 
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with a wave vector k z = kp along the magnetic field returns to the contact when the z- 
component of the momentum k zs = ZQm*Q/2irsh, for integer s. For these orbits the time 
of the motion over a distance zq in the z direction is a multiple of the cyclotron period 
Tjj = 2tt/Q. Thus, after s revolutions the electron returns to the contact axis at the point 
z = 0. The phase which the electron acquires along the spiral trajectory is composed of two 
parts, A0 = A0-L + A0 2 . The first, A(p 1 = k zs z is the 'geometric' phase accumulated by an 
electron with wave vector k zs over the distance z . The second, A0 2 = ns(eHr 2 /ch) is the 
phase acquired during s rotations in the field H, where r s = hc\/k\ — k 2 s /eH is the radius 
of the spiral trajectory. Substituting k zs and r s in the equation for A0 we find 

A0 = 2ns£ F /hn + z 2 /4iisa 2 H . (70) 

This is just the phase shift that defines the period of oscillation in the contribution AG 2 
( 169]) to the conductance. It describes a trajectory which is straight for the part from the 
contact to the defect and spirals back to the contact by s windings as it is shown in Fig. 8. 
There are trajectories consisting of helices in the forward and reverse paths, with s and 
s' coils, respectively. However, the contribution of these trajectories to the conductance is 
smaller than AG 2 (16 9 p by a factor ~ 1/ (k-pan) "C 1. Note that, although the amplitude of 
the oscillation AG 2 (169 j) is smaller by a factor Ml/e-p than the amplitude of the contribution 
AGi ([HSD, the first depends on the depth of the defect as Zq 3 ^ 2 while AGi ~ z^ 2 . The slower 
decreasing of the amplitude for AG2 is explained by the effect of focusing of the electrons in 
the magnetic field. The predicted oscillations, Eq. ( )69|) . are not periodic in H nor in 1/H. 
Their typical period can be estimated as the difference AH between two nearest-neighbor 
maxima 

The period (1711) depends on the position of the defect. It is larger than the period of de 
Haas- van Alphen oscillation, (AH/H) dHvA ~ hQ/ep. Both of these periods are of the same 
order of magnitude. 



VII. NONMAGNETIC DEFECT IN A SUPERCONDUCTOR 

In this section we present the results of a theoretical investigation of the conductance, 
G ns , of a normal metal - superconductor (NS) point contact (with radius a < \p ) in the 
tunnelling limit and discuss the quantum interference effects originating from the scattering 



of quasiparticles by a point-like nonmagnetic defect j47j. The model is described in Sec. II 
and illustrated in Fig. [H modified with having the half-space z > occupied by a (s-wave) 
superconductor. At zero temperature a tunnel current flows through the contact for an 
applied bias eV larger than the energy gap of the superconductor A . In order to evaluate 
the total current through the contact, I (V) , the current density jk (r) of quasiparticles 
with momentum k at z > 0, formed by electrons transmitted through the contact must 
be found. The current density jk (r) is expressed in terms of the coefficients (r) and 
Vk (r) of the canonical Bogoliubov transformation | 70| . The functions (r) and (r) 
satisfy the Bogoliubov-de Gennes (BdG) equations [71] . which must be supplemented with 
a self-consistency condition for the order parameter A (r), and boundary conditions which 
connect «k and t>k in the normal metal to those in the superconductor at the contact. For 
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a tunnel contact one can neglect Andreev reflections, because these lead to corrections to 
the conductance proportional to |i| 4 (72[; the functions «k and t>k satisfy the same boundary 
conditions (TTOj) . ( TTTj) as the wave function for a contact between normal metals. 

It is obvious that the method described in Sec. II and Sec. Ill can be generalized to NS 
contacts. As a first step the Bogoliubov-de Gennes (BdG) equations must be solved in linear 
approximation in the transmission amplitude t in the absence of the defect, D = 0, after 



an ana- 



which the corrections due to the scattering by the defect can be found. In Ref. [47 
lytical solution for the BdG equations was found with the approximation of a homogeneous 
order parameter A (r) = A 6 (z) . 

At small applied bias eV <C Hud <C £f (wd is the Debye frequency), and in linear 
approximation in the electron-defect interaction constant g the conductance G ns of a NS 
tunnel point contact can be presented as the sum of two terms, 

G (V, r„) = Gons (V) + AG osc (V, r ) , eV > A . (72) 

The first term, Go ns (V), in Eq. f J72|) is the conductance of the NS tunnel point contact in 
the absence of the defect 

G 0ns (V) = G N S (eV) , (73) 

where Go is the conductance of a contact between normal metals (12"71 ). which is multiplied by 
the normalized density of states of the superconductor N s (E) — Ej \/ E 2 — A§ at E = eV 



Although such result is not unexpected and has been confirmed by experiment [22j, for a 



contact of radius k pa < 1 it was not obvious and it is first obtained in Ref. |47J • The second 
term describes the oscillatory dependence of the conductance with the distance between the 
contact and the defect. If the defect is situated in the superconductor (zq < 0) 

AG osc (V jro ) = -G 0ns (V)g(^j J2^a(eV)w(k a r ), (74) 



where 



hi 



V2 



?7r 



h 



e F ±J(eVy-At 



1/2 

2 _ A2 



(76) 



and the function w (k a r ) is given by Eq. ( 1291) . Equation (1741) is obtained by neglecting all 
small terms of the order of Aq/ef and eV/e-F- Note that we kept the second term in square 
brackets in the formula for k± (Eq. f J76|) ) because for large r , (a/ (eV) 2 — Aq/£f)(^f^o) — 1) 
the phase shift of the oscillations may be important. 



VIII. CONCLUSIONS 

Thus, we have reviewed some theoretical aspects of the possibility of investigating sub- 
surface defects by STM experiments. The theoretical results show that the amplitude of the 
oscillations of the STM conductance resulting from quantum interference of electron waves 
injected by the STM tip and scattered by the defect remains sufficiently large (~ 10~ 3 Go)> 
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even for defects located more than 10 atomic layers below the surface. For example, in 
the STM experiments of Ref. [73| signal-to-noise ratios of 5 • 10~ 4 (at 1 nA, 400 Hz sample 
frequency) have been achieved. Recently, the possibility of observing defects at such depths 
below the surface has been demonstrated in experiment j33j . 

The STM tip plays the role of "locator", which detects a defect below the metal sur- 
face by using electron waves. The defect in turn produces information about its (defect) 
characteristics, as well as showing properties of the host metal by producing Friedel-like 
oscillations in the STM conductance. The phase of the oscillations, 2kpTQ, is defined by the 
Fermi wave vector kp and the tip-defect distance vq. One of the possibilities to determine 
the defect depth zq below surface is by changing the maximal value of the wave vector by 
accelerating the electrons with an applied bias eV [5l|, [74| . When the tip is situated above 
the defect, the period of the oscillations in G (V), Ak F (eV) z = ir, uniquely defines z . As 
the period of the oscillations becomes longer for small z the minimum detectable depth will 
be determined by the maximum voltage that can be applied over the junction. For example, 
30 mV is sufficient for probing a quarter of a conductance oscillation caused by a defect at 
1 nm depth. 

Another factor in setting the oscillation phase is the shape of the Fermi surface (FS). 
As was shown, for an anisotropic FS e (k) = Ep the phase and amplitude of conductance 
oscillations depend on the characteristics of the FS in the point for which the direction of 



the velocity v = vo is parallel to the vector ro directed from the STM tip to the defect [55 
Namely, the phase of the oscillations is defined by the projection of k on the direction of 
v , and the oscillation amplitude depends on the curvature of the FS. Depending on the 
geometry of the FS there can be several points with the same direction of the velocity, or, 
if the FS has open parts, certain directions of the velocity can be forbidden. It follows 
from the results above that curves of constant phase kvoro/ |vo| (maxima and minima) in 
the interference pattern of the STM conductance show the contours formed by projections 
of the vector k on the vector normal to the FS. Although such contours reflect the main 
features of the FS geometry, they cannot be considered as a direct imaging of the FS. 

Electron scattering by subsurface magnetic defects in STM conductance possesses some 
features distinct from the scattering by magnetic adatoms and the form and sign of the 
Kondo anomaly due to a subsurface magnetic defect depend on the depth ^|. Near the 
Kondo resonance the scattering phase shift Sq tends to tt/2, and including multiple electron 
scattering events after reflections by the metal surface becomes essential. This explains the 
appearance of harmonics in the oscillatory part of the conductance, which have an additional 
phase shift A0 = 2 (n — 1) kpzo + nSo, where n is the number of electron reflections by the 
surface. The determination of this phase shift near the Kondo resonance (V ~ Vk) and far 
from it (where 5q <C 1,) for the first (n = 0) and second (n = 1) harmonics provides an 
alternative way to find the depth of the defect z . 

The possibilities of investigating magnetic defects are extended by injecting a spin- 
polarized current. If the subsurface cluster possesses an unscreened magnetic moment ^ cff , 
the scattering amplitudes of spin-up and spin-down electrons are different. This results 
in a dependence of the oscillation amplitude on the angle between the vector /j, eS and the 
polarization direction of the STM current - referred to as a magneto-orientational effect 67 . 

A strong magnetic field perpendicular to the metal surface changes the interference pat- 
tern in the STM conductance fundamentally. As a result of Zeeman splitting ±g e fi B H of the 
Landau energy levels the interference patterns formed in the dependence G (ro) by electrons 
with different spin directions do not coincide because of different electron wave lengths for 
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energies Ef ± g e /j, B H. The superposition of the two oscillatory parts may result in a beating 
of the total amplitude of the oscillations. Along with the well-known quantum oscillations 
having the periodicity in H~ l of the de Haas - van Alphen effect in the STM conductance, 
in the presence of a defect two new types of oscillations are present. The first is related to 
flux quantization through the projection of the electron trajectory on the surface plane. The 
second type of oscillation G (H) is related to a focusing effect of the magnetic field. As in 
Sharvin's two-point contact experiments, in which electrons were focused on a collector by 
a magnetic field directed along the line connecting the contacts (geometry of longitudinal 
electron focusing) j75|, the magnetic field can periodically focus the electrons injected by 
the tip onto the defect. This results in periodic increasing or decreasing of the part of the 
conductance related to the scattering by the defect (66| . 




If the electrons tunnel from a normal-metal STM tip into a superconductor the wave 
incident on the contact is transformed into a superposition of 'electron-like' and 'hole-like' 
quasiparticles. In the case of a location of the defect in the superconductor quantum in- 
terference takes place between the partial wave that is transmitted and the one that is 
scattered by the defect, for both types of quasiparticles independently (Eq. (]74"|)). Although 
the difference between wave vectors k^' (eV) of 'electrons' and 'holes' is small the shift 
— k^) r between the two oscillations should be observable |47j ]. 
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